Single Cell Molecular Heterogeneity In Musculoskeletal Differentiation

Mesenchymal stem cells (MSCs) display substantial cell-to-cell variation that manifests across many aspects of cell phenotype and complicates the use of MSCs in regenerative applications. However, most conventional assays measure MSC properties in bulk and, as a consequence, mask this cell-to-cell variation. To better understand MSC heterogeneity and its underlying mechanisms, we quantitatively assessed MSC phenotype within the context of chondrogenesis, amongst clonal populations and single cells. Clonal MSCs differed in their contractility, ability to transmit extracellular strain the nucleus, capacity to form cartilage-like matrix, and transcriptomic signature. RNA FISH measurements of single cell gene expression found that both primary chondrocytes and chondrogenically-induced MSCs showed substantial mRNA expression heterogeneity. Surprisingly, variation in differentiation marker transcript levels only weakly associated with cartilage-like matrix production at the single cell level. This finding suggested that, although canonical markers have very clear functional roles in differentiation and matrix formation, their instantaneous mRNA abundance is only tenuously linked to the chondrogenic phenotype and matrix accumulation at the single cell level. One possible explanation for the apparent disconnect between gene and protein expression is that mRNA and protein exhibit different temporal dynamics. Using stochastic models of single cell behavior, we explored the impact of transcriptional stochasticity and temporal matrix dynamics on the perceived relationship between single cell mRNA and protein abundance. Simulations suggested that considering recent temporal fractions of protein (vs. total protein) increased the correlation between mRNA and protein abundance, and illustrated that mRNA stability was a crucial determinant of the timescale over which any such correlation persisted. Experimentally, non-canonical amino acid tagging was used to visualize and quantify temporal fractions of nascent extracellular matrix with high fidelity. The organization and temporal dynamics of the proteinaceous matrix depended on the biophysical features of the microenvironment, including the biomaterial scaffold and the niche constructed by cells themselves. Both chondrocytes and MSCs demonstrated marked cell-to-cell heterogeneity in nascent matrix production, consistent with model predictions. Ongoing work aims to combine these experimental measurements of nascent protein expression with instantaneous measures of mRNA abundance to better understand the mRNA-protein relationship, and to harness this understanding to improve regenerative therapies

Genetic Markers of Toxicity From Capecitabine and Other Fluorouracil-Based Regimens: Investigation in the QUASAR2 Study, Systematic Review, and Meta-Analysis

Fluourouracil (FU) is a mainstay of chemotherapy, although toxicities are common. Genetic biomarkers have been used to predict these adverse events, but their utility is uncertain

Identification of unique neoantigen qualities in long-term survivors of pancreatic cancer

Pancreatic ductal adenocarcinoma is a lethal cancer with fewer than 7% of patients surviving past 5 years. T-cell immunity has been linked to the exceptional outcome of the few long-term survivors1,2, yet the relevant antigens remain unknown. Here we use genetic, immunohistochemical and transcriptional immunoprofiling, computational biophysics, and functional assays to identify T-cell antigens in long-term survivors of pancreatic cancer. Using whole-exome sequencing and in silico neoantigen prediction, we found that tumours with both the highest neoantigen number and the most abundant CD8+ T-cell infiltrates, but neither alone, stratified patients with the longest survival. Investigating the specific neoantigen qualities promoting T-cell activation in long-term survivors, we discovered that these individuals were enriched in neoantigen qualities defined by a fitness model, and neoantigens in the tumour antigen MUC16 (also known as CA125). A neoantigen quality fitness model conferring greater immunogenicity to neoantigens with differential presentation and homology to infectious disease-derived peptides identified long-term survivors in two independent datasets, whereas a neoantigen quantity model ascribing greater immunogenicity to increasing neoantigen number alone did not. We detected intratumoural and lasting circulating T-cell reactivity to both high-quality and MUC16 neoantigens in long-term survivors of pancreatic cancer, including clones with specificity to both high-quality neoantigens and predicted cross-reactive microbial epitopes, consistent with neoantigen molecular mimicry. Notably, we observed selective loss of high-quality and MUC16 neoantigenic clones on metastatic progression, suggesting neoantigen immunoediting. Our results identify neoantigens with unique qualities as T-cell targets in pancreatic ductal adenocarcinoma. More broadly, we identify neoantigen quality as a biomarker for immunogenic tumours that may guide the application of immunotherapies

An original phylogenetic approach identified mitochondrial haplogroup T1a1 as inversely associated with breast cancer risk in BRCA2 mutation carriers

Introduction: Individuals carrying pathogenic mutations in the BRCA1 and BRCA2 genes have a high lifetime risk of breast cancer. BRCA1 and BRCA2 are involved in DNA double-strand break repair, DNA alterations that can be caused by exposure to reactive oxygen species, a main source of which are mitochondria. Mitochondrial genome variations affect electron transport chain efficiency and reactive oxygen species production. Individuals with different mitochondrial haplogroups differ in their metabolism and sensitivity to oxidative stress. Variability in mitochondrial genetic background can alter reactive oxygen species production, leading to cancer risk. In the present study, we tested the hypothesis that mitochondrial haplogroups modify breast cancer risk in BRCA1/2 mutation carriers. Methods: We genotyped 22,214 (11,421 affected, 10,793 unaffected) mutation carriers belonging to the Consortium of Investigators of Modifiers of BRCA1/2 for 129 mitochondrial polymorphisms using the iCOGS array. Haplogroup inference and association detection were performed using a phylogenetic approach. ALTree was applied to explore the reference mitochondrial evolutionary tree and detect subclades enriched in affected or unaffected individuals. Results: We discovered that subclade T1a1 was depleted in affected BRCA2 mutation carriers compared with the rest of clade T (hazard ratio (HR) = 0.55; 95% confidence interval (CI), 0.34 to 0.88; P = 0.01). Compared with the most frequent haplogroup in the general population (that is, H and T clades), the T1a1 haplogroup has a HR of 0.62 (95% CI, 0.40 to 0.95; P = 0.03). We also identified three potential susceptibility loci, including G13708A/rs28359178, which has demonstrated an inverse association with familial breast cancer risk. Conclusions: This study illustrates how original approaches such as the phylogeny-based method we used can empower classical molecular epidemiological studies aimed at identifying association or risk modification effects.Peer reviewe

Finishing the euchromatic sequence of the human genome

The sequence of the human genome encodes the genetic instructions for human physiology, as well as rich information about human evolution. In 2001, the International Human Genome Sequencing Consortium reported a draft sequence of the euchromatic portion of the human genome. Since then, the international collaboration has worked to convert this draft into a genome sequence with high accuracy and nearly complete coverage. Here, we report the result of this finishing process. The current genome sequence (Build 35) contains 2.85 billion nucleotides interrupted by only 341 gaps. It covers ∼99% of the euchromatic genome and is accurate to an error rate of ∼1 event per 100,000 bases. Many of the remaining euchromatic gaps are associated with segmental duplications and will require focused work with new methods. The near-complete sequence, the first for a vertebrate, greatly improves the precision of biological analyses of the human genome including studies of gene number, birth and death. Notably, the human enome seems to encode only 20,000-25,000 protein-coding genes. The genome sequence reported here should serve as a firm foundation for biomedical research in the decades ahead

DNA methylation patterns identify subgroups of pancreatic neuroendocrine tumors with clinical association

Here we report the DNA methylation profile of 84 sporadic pancreatic neuroendocrine tumors (PanNETs) with associated clinical and genomic information. We identified three subgroups of PanNETs, termed T1, T2 and T3, with distinct patterns of methylation. The T1 subgroup was enriched for functional tumors and ATRX, DAXX and MEN1 wild-type genotypes. The T2 subgroup contained tumors with mutations in ATRX, DAXX and MEN1 and recurrent patterns of chromosomal losses in half of the genome with no association between regions with recurrent loss and methylation levels. T2 tumors were larger and had lower methylation in the MGMT gene body, which showed positive correlation with gene expression. The T3 subgroup harboured mutations in MEN1 with recurrent loss of chromosome 11, was enriched for grade G1 tumors and showed histological parameters associated with better prognosis. Our results suggest a role for methylation in both driving tumorigenesis and potentially stratifying prognosis in PanNETs